General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 








NASA CH- 167926 
BAT Report #05636-941009 


NASTRAN FLUTTER ANALYSIS 
OF 

ADVANCED TURBOPROPELLERS 


by 

V. Elchuri 
G. C, C, Smith 


BELL AEROSPACE TEXTRON 
P, 0. Box One 
Buffalo, New York 14240 


iNASl-CB«IB7926) NAS2BAB PIBSOS AifU.X3IS K84- 14148 

0^ AnYAHCBP rcUBBOPECEBXIEIiS EinaX leport 
(Textron Bell Aerospace Co,, Buffalo, N, X,) 

121 p HC A06/MP A01 CSCl, 02A tinclas 

G3/07 42838 



National Aeronautics and Space Administration 
Contract NAS3-22533 


NASA Lewis Research Center 
Clevelands Ohio 44135 


"April 1982 



\. Heport No. 2 Gowr.mcni Acctfsipn No. 

_ fmSA CR-167926 

3. ficcipiCnl’t Ol.Jlog No 

*1. Titfc and Subtilif* 

NASTRAN MUTTER AMAlYEIS OF ADVANCED TURBOPROPELLERS (U) 

S, (U’pori Data 

April 1982 

G. Pcrformmg Orgiitualion Code 

7. Authofisl 

V. Elchuri and G. C. C. Smith 

0. Porforr.iing OtganitQtion Pepoft No. 
02y36-g<U009 

to Work Unit No 

D. Pcf forming Or.'VUiatiDn Nama and Addreu 

Bell Aerospace Textron 
P. 0. Box One 
Buffalo, New York imo 

11. Coniract or Grant Nts 
NAS 3-22533 

13. Type of nojiort and Period Covered 
Contractor Report 

12. SponsorinQ Agency Nanio and AUdruis 

National Aeronautics and Space Administration 
Nasliington, D. C, 20bf1G 

M. Spantaring Ageney Code 

1 

15. Sup;ilcmcniofv Nnles 

NASA Program Nanagrr - Richard E. Morris » Final Report 

10 Abslrncl 


An G listing capability developod by the authors to conduct modal flutter analysis of tuned 
bladad-shrouded discs in NASTP.AN has been m''‘’"iod and applied to investigate the subsonic 
unstalled flutter characteristics of advanced t.prbopropellers. 


The modifications pertain to the inclusijn of oscillatory modal aerodynamic loads of blades 
with large (backward and forward) variable swetp. The two-dimensional subsonic cascade 
unsteady aerodynamic theory of W. P, Jones and 0, M. Rao has been applied in a strip theory 
manner with appropriate modifications for tl>- sweep effects. Each strip is associated with 
a chord selected normal to any spanwise reference curve such as the blade leading edge. 

The stability of three operating conditions of a 10-bladed propeller is analyzed. Each of 
these operating conditions is iterated once to determine the flutter boundary. A S-bladed 
propeller is also analyzed at one operating condition to investigate stability. Analytical 
results obtained are in very good agreement with those from wind tunnel tests. 

The work was conducted under Contract MAS 3-22532 from NASA Lewis Research Center, Cleveland, 
Ohio, with Hr. Richard E. Morris as the Technical Monitor. 


17. Key V7ords (Suggested by Aullior{s)l 

Advanced Turbopropellers, NASTRAN, Flutter 
Analysis, Finite Elements 

IB. Distribution Statement 
Unlimited - Unclassified 
Subject Categories: 07,61 

19, Security Ciiutlf. lot tlili report) 

Unclassified - Unlimited 

20. Security Classif, (of IhlJ pagcl 
Unclassified - Unlimited 

2t. No. of PngoE 
112 

22. Price* 


‘ Tor sale by the Nalloiial Technical Iniofiiialion Setvice, Spiinelield, Vitginia 2ZlSi 


NASA-r-lfiB (Kev. 10-7S) 


i 



NASA CR- 167926 
BAT Report ^02536-941009 


NASTRAH FLUTTER ANALYSIS 
OF 

ADVANCED TURBOPROPELLERS 


by 

V. Elchuri 
G. C. C. Smith 


BELL AEROSPACE TEXTRON 
P. 0. Box One 
Buffalo, New York 14240 


National Aeronautics and Space Administration 
Contract NAS3-22533 

NASA Lewis Research Center 
Cleveland, Ohio 44135 


April 1982 


ABSTRACT 


An existing capability developed by the authors to conduct modal flutter 
analysis of tuned bladed-shrouded discs in NASTRAN has been modified and applied 
to investigate the subsonic unstalled flutter characteristics of advanced 
turbopropellers. 

The modifications pertain to the inclusion of oscillatory modal aerodynamic 
loads of multi-bladed propellers with large (backward and forward) variable 
sweep. Kaza and Kielb of NASA Lewis have demonstrated the aerodynamic signifi- 
cance of cascade and sweep effects on simplified structural models. In the 
current work, the two-dimensional subsonic cascade unsteady aerodynamic theory 
of W. P. Jones and B. M. Rao has been applied in a strip theory manner with 
appropriate modifications for the sweep effects. Each strip is associated with 
a chord selected normal to any spanwise reference curve such as the blade leading 
edge. 

The sUbility of three operating conditions of a 10-bladed propeller is 
analyzed. Each of these operating conditions is iterated once to determine the 
flutter boundary. A 5-bladed propeller is also analyzed at one operating condi- 
tion to investigate stability. Analytical results obtained are in very good 
agreement with those from wind tunnel tests. 

The work was conducted under Contract NAS3-22533 from NASA Lewis Research 
Center, Cleveland, Ohio, with Mr. Richard E. Morris as the Technical Monitor. 


SUMMARY 


Generalized aerodynamic forces for varying sweep turbopropeller blades 
have been developed based on subsonic two-dimensional cascade theory. A 
bladed-disc dynamics/flutter program has been modified to enable the study 
of adviinced turbopropeller flutter. Specific analyses show excellent agree- 
ment with wind tunnel results. 
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1.0 INTRODUCTION 


As part of the government/ Indus try effort in evolving aircraft of high 
efficiency, NASA Lewis has been involved in the study ot advanced turbopro- 
pellers. These devices have many thin blades of a complex shape as compared 
with conventional propellers. As is historically frequent with advanc^^d 
aeronautical lifting surfaces, the Advanced turbopropeller has shown evidence 
of flutter problems in wind tunnel tests. In order to compliment and support 
wind tunnel evidence and other analytical studies, Bell Aerospace Textron has, 
under NASA sponsorship, modified a computer program previously developed for 
Bladed Disc Dynamics and applied it to studies of Advanced Turbopropeller 
Flutter. 

This report describes the theoretical background for those modifications 
and the specific applications and results obtained. The details of the modi- 
fications necessary for generalized aerodynamics of varying sweep blades are 
given in Appendices. 

The necessary coding modifications to the aerodynamic modules of the 
existing program have been documented and the program and documentation have 
been delivered to NASA Lewis. The program is operational in the NASTRAN 
General Purpose Computer Program at Level 17.7 on IBM and UNIVAC systems. 
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2.0 THEORY 


Mu1t1-b1aded advanced turbopropellers are geometrically cyclic structures with 
low aspect ratio thin blades of varying sweep. The blades are mounted on a 
relatively rigid hub (Figure 1) and, therefore, can be considered to be struc- 
turally independent. This permits modal analysis of only one root-fixed blade 
without recourse to special harmonic analysis techniques applicable to cyclic 
structures. From a flutter aerodynamics viewpoint, tne estimation of the 
generalized oscillatory aerodynamic loads on the propeller blades depends on 
the aerodynamic theory employed. In the present work, the two-dimensional sub- 
sonic cascade unsteady aerodynamic theory of Jones and Rao (Ref. 1) is applied 
in a strip theory manner with appropriate modifications recognizing the vari- 
ability of blade sweep with radius. 

The following sections di-it iss the structural and aerodynamic modelling, 
the equations and the method of modal flutter analysis of advanced turbopropellers. 

NASTRAN Structural and Aerodynamic Models 

With the assumption of a rigid hub and identical (tuned) blades, it suffices 
to model only one blade as shown in Figure 2. The finite element capabilities 
for structural modelling provided in the NASTRAN general purpose finite element 
program are used. 

To facilitate the use of a two-dimensional cascade unsteady aerodynamic 
theory, the aerodynamic model is based on a grid defined by the intersection of 
a series of chords and "computing stations" as shown by the thick solid lines 
in Figure 2. 

In order to apply strip theory in a manner similar to that of Ref. 2, 
the chords are selected normal to any "spanwise" reference curve such as the blade 
leading edge. The choice of the number and location of the chords and the 
computing stations is dictated by the expected variation of the relative flow 
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Figure 1 . An Advanced Turoopropeller 
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Figure 2. NASTRAN Structural and Aerodynamic Models of 

the Advanced Turbopropeller for Flutter Analysis 
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properties across the blade span, and the complexity of the mode shapes exhibited 
by the p* .peller blade. Due to its resemblance to the structural model of the 
blade, and the adequacy of a relatively coarse grid to describe the spanwise 
flow variations, the aerodynamic model is chosen as a subset of the structural 
model as shown in Figure 2. 


Equations of Motion for Modal Flutter Analysis 

If {u} represents the physical degrees of freedom of any one representative 
blade, the equations of motion for flutter analysis of the advanced turbopropeller 
can be written as 


[Ml (U] + (1 + ig) 


[K®] + m 


{u} - [Ql {u} = {0}, 


( 1 ) 


where M is the mass matrix, K® is the elastic stiffness matrix, K'^ and Q are 
the differential stiffness matrix and the aerodynamic coefficient matrix, 
respectively, at the operating condition being considered. For the KE method 
of flutter analysis (Ref. 3 ) used in the present work, an artificial structural 
damping parameter g is introduced in the flutter equations. Although the assump- 
tion of harmonic motion in deriving the aerodynamic coefficients limits the 
flutter equations and their solutions to be strictly valid only at g=0, g->0- 
is taken to represent a stable system approaching flutter. 

For modal flutter analysis, equations (1) are reformulated with the undamped 
structural modes [(pj, computed with the effects of steady state loads taken into 
account, as the bases: 

+ (1 + ig)[K,,] - [Q,,]] (X,) = (0). 


(2) 


Where 


(u) = {^} e’“* = (x,) 


= wfww. 

[K.,] = [♦f[[K‘’] + [k‘'] 

= [♦]’'[Q]W . 


[4>], and 


(3) 
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A detailed derivation of the generalized aerodynamic coefficients matrix 
for the swept blades of advanced turbopropellers Is presented In Appendix A. 
For the KE method of flutter analysis In NASTRAN, equations (2) are rewritten 
as a complex elgenvalu'^ problem 

[cm,,] ♦ [K, ,]]({,)= {OK (4) 


where 


p = - V- (f + 1) Is the complex eigenvalue, 

/l + 1g 


° CQ,,] . and 

^,ref ^s.ref 


k- 


.ref 


V-ref 


Is the reduced frequency. 


^ (5) 


Method of Flutter Analysis 

An overall flowchart for modal flutter analysis of advanced turbopropellers 
is shown In Figure 3. 

For an N-bladed propeller with the blades set at a known reference pitch 
angle, the operating condition to be examined for stability is defined in terms 
of rotational speed and free stream density, velocity and Mach number. 

The steady state loads are used to estimate the differential stiffness of 
the propeller blade. In the examples presented In this report, the differential 
stiffness due to only centrifugal loads is estimated by iteratively solving 

[Ck'] + CK‘'(u)]j{u) . p'f (6) 

for I. to determine K^(u). 

The real eigenvalue problem 


-..^[M] t [CK=] + [K=‘] 
- 


{U} = {0} 


(7) 
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Figure 3. Overall Flowchart of Advanced Turbopropeller Modal Flutter 
Analysis (continued). 
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Figure 3. 


Overall Flowchart of Advanced Turbopropeller 
Modal Flutter Analysis (Concluded) 








Is solved to determine the natural frequencies and mode shapes of the propeller 
blade. 

For selected values of reduced frequency, k- and all possible values 
of the inter-blade phase angle, o, a series of generalized oscillatory aero- 
dynamic coefficients matrices is computed as shown in Appendix A. 

For each of the possible combinations of (p^ k- ^g^, o), the complex 
eigenvalues of the flutter equations, (4), are determined. The complex eigen- 
value 


yields 

“V.ref * • 

g = 2 Re(p)/V- ^g^, and 


y ( 8 ) 


r _ 1 *^s,ref ’ ^s,ref 

2^' W ■ J 

For each {pj ^g^, o) combination, Vj ^g^-g and V- ^g^-f plots of the type 
shown in Figure 4 are generated. At the {p^ ^g^, V- ^g^) representing the 
operating condition being examined, the stability of each of the structural 
modes included for flutter analysis is examined by the value of the damping 
parameter g being ^0. g>0 indicates an unstable mode. In Figure 4, mode 1 

(first bending) with an in-vacuo frequency of 187.9 Hz. (V- ^g^ = 0) is seen 
to be unstable at the operating condition V- ^g^ = 7046 in/sec. at an inter- 
blade phase angle of -108 degrees. 

The root locus of this mode with inter-blade phase angle as the parameter, 
at the (A- V-r representing the operating condition is shown in Figure 5. 
The ordinate is the non-dimensional frequency v = f/f„^^,„,^, and the abscissa is 
the non-dimensional damping p = (g/2)*v. It is seen that this mode is unstable 
for 0 = 72°, 180°, -144°, -108° and -72°, with -108° and -72° representing the 
most and the least unstable conditions, respectively. 
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with regards to the use of two-dimensional subsonic cascade unsteady aero- 
dynamic theories, It Is relevant to mention that for certain combinations of 
Mach number. Inter-blade phase angle, blade spacing to chord ratio and stagger 
angle, there exist critical values of reduced frequency about which the aero- 
dynamic coefficients are very sensitive, and at which they vanish (Ref. 4). A 
brief description of such "acoustic resonance" conditions Is given in Ref. 5. 

In the case of the advanced turbopropeller application, however, the extent to 
which such conditions on isolated chords influence the generalized aerodynamic 
matrix and therefore affect the flutter results. Is difficult to judge. Although 
no attempt has been made In this report to study the influence of the presence 
of acoustically resonant chords, the following points are offered: 

1. At any given operating condition, the ratio of the work done by the aero- 
dynamic forces on the strips associated with all "acoustically resonant" 
chords to the work done by the aerodynamic forces on the entire blade in 
any given structural mode may be a good Indication of the Influence of 
the acoustic resonance. 

2. By increasing the number of chords spanning the blade, the aerodynamic 
work done per strip can be reduced. 

3. Due to the large spanwise variation in the parameters defining the critical 
reduced frequencies, 1t may be possible to eliminate acoustic resonance 
conditions by relocating some of the blade chords. 

The method of flutter analysis discussed in this section can be repeated 
at other operating conditions till neutrally stable structural modes are found. 
The subsonic classical flutter boundary may then be defined. 
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3.0 APPLICATION 


The theoretical development of the preceding section, as Implemented In 
the NASTRAN general purpose finite element program (Ref. 6), was applied to 
conduct the modal flutter analysis of the SR-5 advanced turbopropel ler developed 
by the NASA Lewis Research Center and the Hamilton Standard Division of the 
United Technologies Corporation. Two (five- and ten-bladed) propellers were 
analyzed. 

The four operating conditions for analysis were selected by NASA LeRC from 
wind tunnel test conditions, and are shown In Table 1. 

The propeller hub was considered rigid compared to blades, and only one 
of the blades was modelled. The NASTRAN structural model, shown In Figure 2, 
consisted of 130 bending-cum-membrane plate elements with a total of 155 grid 
points. 

The aerodynamic model, also shown in Figure 2, consisted of 10 chords with 
four computing stations each. The chords were selected approximately normal to 
the blade leading edge. For the four operating conditions, the corresponding 
unsteady aerodynamics data are presented In Tables 2 through 5 . For case 3, 
the rate of change of blade chord with spanwise reference distance, 9cj/9z, was 
computed from a finer (20 chords) mesh. 

To investigate flutter instability, the following steps were taken at 
each of the operating conditions; 

1. Centrifugal loads at the operating rpm were applied to the as manufactured 
(pretwisted) blade shape to compute the differential stiffness. 

2. The displacements at the end of the differential stiffness calculations 
were added to the pretwisted blade shape to define the deformed blade 
shape. The elastic and differential stiffnesses were then used to deter- 
mine the natural frequencies and mode shapes of the deformed blade. 
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TABLE 3. UNSTEADY AERODYNflMICS DATA FOR CASE 2 
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TABLE 5. UNSTEADY AERODYNAMICS DATA FOR CASE 4 





















































3. The first six structural modes of the deformed blade were used to derive 
the generalized aerodynamic forces on the blade. The KE method of modal 
flutter analysis was used to obtain the subcritical, critical (flutter) and 
supercritical eigenvalues and to generate the V>g and V-f curves. 

4. The V-g curves were interpreted to determine the stability of the structural 
modes. For the unstable modes, the root locus plots with inter-blade phase 
angle as the parameter were analyzed to determine the dominant flutter 
inter-blade phase angle. 

For the 10-bladed propeller (Table 1 , cases 2-4), steps 3 and 4 were 
repeated once in order to locate the neutrally stable operating conditions 
(flutter boundary). In order to retain the structural modes computed in step 2, 
rpm's were unchanged while varying the free stream Mach numbers. 

The results of the analysis are presented and discussed in Section 4. 


4.0 RESULTS AND DISCUSSION 


Results of differential stiffness, modal and flutter analyses for all the 
cases discussed in the preceding section are presented with a view to 

1. determine if the given operating conditions (cases 1 through 4) are stable 
or unstable, and 

2. locate the subsonic classical flutter boundary for the 10-bladed propeller 
(cases 2 through 4, and their iterations). 

Figure 6 schematically illustrates the streamwise blade tip deflections 
obtained at the end of the differential stiffness computations for cases 1 
through 4. These computations included the effects of only centrifugal loads 
and neglected those due to steady state airloads. The blade tip, represented by 
a straight chord, is seen to deflect 4.3 degrees to 5.8 degrees as the propeller 
speed varies from 6000 to 6900 rpm. 

Figures 7 and 8 are for case 3, and f/pically show the deformed blade 
shapes obtained by excluding and including the differential stiffness, respec- 
tively. 

Figure 9 presents a summary of the modal analysis of the SR-5 blade for 
cases 1 through 4 and static condition (0 rpm). The differential stiffness was 
included for each case in determining the natural frequencies and mode shapes. 
Bench test frequencies available for the first four natural modes, under static 
conditions, are also shown in the figure. Excellent agreement with analysis is 
indicated except for the third mode frequency which is predicted about 16% 
higher. As seen in Figure 12 , this is primarily the first torsion mode about 
the blade pitch axis, and the higher frequency is possibly due to the over- 
representation of the torsional stiffness by using plate elements near the 
blade root and shank. 
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SR-5 STREAMWISE Bl ADE TIP DEFLECTIONS 



Figure 7. SR-5 Steady State Deflections Without Differential Stiffness 




Figure 8. SR-5 Steady State Deflections with Differential Stiffness 

(Case 3} 
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Figures 10 through 15 , and 16 through 21 , respectively, represent 
the SR-5 blade mode shapes under static and case 3 conditions. Table 6 lists 
the contour values for these figures. Figures 10 through 13 also show the 
mode shapes determined by holographic bench tests. The comparison between the 
computed and observed mode shapes is generally good. 

Figures 22 through 26 (case 1), 27 through 36 (case 2), 37 through 46 
(case 3} and 47 through 56 (case 4) are the V-g, V-f curves generated by the 
solution of the flutter equations as a complex eigenvalue problem. For each of 
these figures, 

1. the density of the inflow was constant, and as given by Table 1 , 

2. the inter-blade phase angle was constant at one of its permissible values: 

o = 2tt n/N , 

with N = the number of blades on the propeller, 

n = 0, tl, t2 ±Np N 2 » 

where for N odd , 

N-| = (N-l)/2, ^2 does not exist, and 
for N even , 

= (N-2)/2, = N/2, 

3. the reduced frequency at reference chord, k- was varied as 0.1, 0.2, 
0.3, 0.6, 0.9, 1.2 and 1.5 to cover the frequency range of the first six 
structural modes, 

4. a vertical line at the V— correspono ing to the operating condition is 
drawn to facilitate determination of the equilibrium of the various 
structural modes from the V-g curves. 

With regards to the V-g, V-f curves for all four cases, the following 
general observations are made: 

1. At the given operating conditions, the frequencies of modes 1, 3 and 5 

increased, with mode 1 gaining the most. The frequencies of modes 2, 4 and 
6 decreased marginally. 


25 



ORIGINAL PAQ* W 
OF POOR QUALITY 


Figure 10 . SR-5 Mode Shapes, 0 rpm. Mode 1 
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Figure 13. SR-5 Mode Shapes, 0 rpm. Mode 4 







Figure 15. SR-5 Mode Shape, 0 rpm, Mode 6 
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Table 6. CONTOUR VALUES FOR MODAL PLOTS 
(Figures 10-21) 


Contour Symbol 

Y-Displacement 
(Figures 10-15) 

Y-Displacement 
(Figures 16-21) 

1 

-0.003 

-1.0 

2 

0.11 

-0.9 

3 

0.22 

-0.78 

4 

0.33 

-0.67 

5 

0.44 

-0.56 

6 

0.55 

-0.44 

7 

0.67 

-0.33 

8 

0.78 

-0.22 

9 

0.89 

-0.10 

10 

1.0 

+0.005 


Notes ; 

1. For Figures 10 through 21 : 

a) Contours are plotted on the undeformed shape. 

b) Y axis is normal to the axis of rotation, and approximately 
normal to the blade mean surface. 

2. For Figures 10 through 15 ; 

a) Outline of the deformed shape is also shown. 

b) Contour with symbol approximately represents a nodal line. 

3. For Figures 16 through 21 : 

a) Contour with sumbol '10' approximately represents a nodal line. 
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Figure 55. V-g, V-f Curves, Case 4 
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2. At the given operating conditions, only mode 1 (first bending mode repre- 
sented by the symbol Q ) shows instability for a limited number of inter- 
blade phase angles. The degree of instability, +g, varies with the inter- 
blade phase angle. 

3. The critical reduced frequency lies between 0.2 and 0.3. 

For each of the four cases at their respective operating conditions, 

Tables 7 through 10 summarize the root locus of the first bending mode with 
inter-blade phase angle as the parameter. A non-dimensional frequency, 

a non-d( t-sional damping, 

\1 s (g/2) • V 

have been defined to plot the root locus shown in Figures 57 through 60 (re- 
presented by the symbol q ) . 

In these figures, though the results are valid only at the discrete Inter- 
blade phase angles identified by symbols, 'smooth' curves have been drawn for 
clarity of presentation. The curve-fitting routine used for plotting these 
curves was supplied with the zero inter-blade phase angle as the starting and 
end points to define the closed loops, and explains the sudden change in curva- 
ture about that point. 

With regards to Figures 57 through 60 , it is observed that; 

1. The number of inter-blade phase angles for which the first bending mode 
is unstable, varies from case to case. 

2. Since the m= 0 line represents neutral equilibrium, the degree of instability 
indicated by the positive values of m can be used to estimate the location 
of the flutter boundary. Of the 10-bladed propeller cases (2, 3 and 4), 
cases 3 and 4 seem to be closest to and farthest from the flutter boundary. 

3. For the 10-bladed propeller, the root behavior at o=72° is more sensitive 
to changes in operating conditions than at other inter-blade phase angles. 
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TABLE 7 


ROOT LOCUS RESULTS FOR FIRST BENDING MODE, CASE 1 


in/sec. 

'vac. ='89.4 HZ 


j" " ' " '■ "" 1 

Inter-blade 
Phase Angle,. 

0 deg. 

Reduced 

Frequency, 

•^s.ref 

— 

Frequency, 
f Hz. 

Damping, 

9 

Non-dim. 
Frequency , 

““''"vac. 

Non-dim. 
Oampina, 
u = (g/2)*v 

0 

.248 

218.1 

-5.143E-.? 

1.152 

-2.961E-2 

72 

.248 

218.5 

-3.202E-2 

1.154 

-1.847E-2 

144 

.244 

215.0 

5.393E-2 

1.135 

3.061E-2 

•144 

.246 

216.4 

6.440E-2 

1.143 

3.679E-2 

-72 

.237 

209.1 

-2.781E-2 

1.104 

-1.535E-2 





























TABLE 8. ROOT LOCUS RESULTS FOR FIRST BENDING MODE, CASE 2 

V- = 7495 in/sec 

f... = 175.9 Hz. 


Inter-blade 
Phase Angle, 
a deg. 


Reduced 

Frequency, 

’^s.ref 


Frequency, 
f Hz. 


Damping 

9 


Non-dim. 

Frequency, 


192.3 


-8.839E-2 


1.093 


Non-dim. 
Damping, 
u ® (g/2)*v 


-4.832E-2 


36 

.235 

192.6 

-8.007E-2 

1.095 

-4.384E-2 

72 

.242 

198.4 

1.698E-2 

1.128 

9.576E-3 

108 

.246 

201.8 

-6.165E-3 

1.147 

-3.536E-3 

144 

.248 

203.2 

1.980E-3 

1.155 

1.144E-3 

180 

1 

.247 

202.6 

1.208E-2 

1.152 

6.957E-3 

-144 

.245 

201.0 

1.823E-2 

1.143 

1 .042E-2 

-108 

.242 

1 

198.6 

1.862E-2 

1.129 

1.051E-2 

-72 

.238 1 

195.8 

1.269E-2 

1.113 

7.063E-3 

-36 

.226 ; 

185.0 

4.207E-3 

1.052 

2.212E-3 



















































TABLE 9. ROOT LOCUS RESULTS 


Inter-blade 
Phase Angle, 
o deg. 


36 


72 


108 


144 


180 


-144 


-108 


-72 


Reduced 

Frequency, 

k- r 
s,ref 


.265 


.267 


.274 


.277 


.277 


.278 


,277 


.273 


,270 


"s.ref = ’ 




Frequenc 
f Hz. 


204.3 


205.8 


210.8 


213.5 


213.9 


214.3 


213.5 


210.6 


207.8 


-36 


.258 


198.9 



■ 

Non-dim. 

Frequency, 

V = f/t 

'vac. 

Non-dim. 
Oamplna, 
u = (g/2)*v 

-8.380E-2 

1.087 

-4.556E-2 

-7.275E-2 

1.095 

-3.984E-2 

7.761E-3 

1.122 

\353E-3 

-1.253E-2 

1.136 

-7.119E-3 

-5.829E-3 

1.138 

-3.318E-3 

4.364E-3 

1.141 

2.489E-3 

8.805E-3 

1.136 

5.002E-3 

l.OOlE-2 

1.121 

5.610E-3 

5.903E-4 

1.106 

3.264E-4 


-9.077E-3 


-1.715E-2 


1.059 





















































TABLE 10. BOOT LOCUS RESULTS FOR FIRST BENDING MODE. CASE 4 

''s.ref “ ”” 
f .... = 179.0 Hz 


Inter-blade 
Phase Angle, 
a deg. 


Reduced 

Frequency, 

*^s,ref 


Frequency, 
f Hz. 


Non-dim. 

Frequency, 


Non-dim. 

Dampina, 

P = (g/2)*v 


0 

.226 

^ .—Id 

196.8 

-8.708E-2 

1.099 

-4.787E-2 

— l-l 1 ■ ■ N I 

36 

.227 

197.9 

-7.649E-2 

1.106 

-4.228E-2 ' 

72 

.231 

200.6 

5.865E-2 

1.121 

3.286E-2 

— ^ 

108 

.238 

207.0 

3.336E-2 

1.156 

1.929E-2 

144 

.238 

207.7 

4.046E-2 

l."60 

2.347E-2 

180 

.239 

208.0 

4.301E-2 

1.162 

2.499E-2 

-144 

.236 

205.4 

5.309E-2 

1.147 

3.046E-2 

-108 

.233 

203.2 

4.962E-2 

1.135 

2.816E-2 

-72 

.230 

200.1 1 

3.800E-2 

1.118 

2.124E-2 

-36 

.217 

189.3 , 

-1.532E-2 

1.058 

-8.101E-3 

















































ilili 



Dampi 









To Illustrate the remark 2 above, Table 11 has been prepared to estimate 
the operating conditions when g (and hence u) is zero by considering similar 
Inlet velocity triangles at the reference chord as shown In figure 61 • 

For the 10~bladed propeller, based on the given operating conditions and 
those Implied at g-0 from Table 11, three new operating conditions were 
selected (Table 12). In order to use the structural modes computed earlier 
for cases 2, 3 and 4, the propeller rpm's were unchanged while changing the 
tunne’* Mach numbers. 

Analysis of the results of these three 1st Iteration cases, conducted 
similar to that discussed earlier In this section for cases 2, 3 and 4, also 
shows only mode 1 to be closest to neutral equilibrium. Results of these 
three 1st Iteration cases are presented In Tables 13 through 15 , Figures 58 
through 60 (represented by the sumbol A )» and Table 11 . 

Figure 62 presents a summary of all the cases considered for the SR-5 
five- and t*'n-bladed propeller flutter analysis. The subsonic classical 
flutter boundary predicted for the ten-bladed propeller using the analysis 
discussed In this report is also shown in this figure. 

Figure 63 is presented to illustrate the excellent agreement between the 
analytical predictions of the present work and the experimental observations 
of NASA LeRC. The analytical predictions of Kaza and Kielb (Ref. 7 ), approxi- 
mating the structural behavior of the SR-5 blades by simplified, constant-chord, 
constant-sweep equivalent beam model with cascade effects 

are also shown in this figure. 
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TABLE 12. SP.-5 OPERATING CONDITIONS FOR 1st ITERATION {lO-BLADEO PROPELLER) 


Case No. 


Blade 

Setting 

Angle 

Free Stream Conditions 




2, 1st iter. 

I 

69“ 

768.3 

1.9534E-3 

0.71 

3, 1st iter. 

6800 

69“ 

733.6 

2.0244E-3 

0.66 

4, 1st iter. 

6200 

69“ 

759.0 

1.7756E-3 

0.70 
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TABLE 13. ROOT LOCUS RESULTS FOR FIRST BENDING MODE, CASE 2, 1st ITERATION 


r 

I 

j 

j 

i 


"I, ref ' 
fvac. “ 


Inter-blade 
Phase Angle, 
0 deg. 

Reduced 
Frequency , 

^,ref 

Frequency, 
f Hz. 


Non-dim. 

Frequency, 

V ® f/f 

'vac. 


0 

.254 

188.4 

-7.769E-2 

1.071 

-4.161E-2 

36 

.255 

188.8 

-6.895E-2 

1.073 

-3.700E-2 

72 

.262 

194.6 

-1.860E-2 

1.106 

-1.029E-2 

108 

.265 

196.3 

-3.214E-2 

1.116 

-1.793E-2 

144 

.266 

196.9 

-2.568E-2 

1.119 

-1.437E-2 

180 

.265 

196-6 

-1.951E-2 

1.117 

-1.090E-2 

-144 

.264 

195.0 

-1.441E-2 

1.113 

-8.020E-3 

-108 

.262 

194.1 

-1.158E-2 

1.103 

-6.389E-3 

-72 

.258 

191.6 

-9.723E-3 

1.089 

-5.295E-3 

-36 

.251 

185.7 

4.066E-3 

1.056 

2.146E-3 





























































TABLE 14. ROOT LOCUS RESULTS FOR FIRST BENDING MODE, CASE 3, 1st ITERATION 


V- ■ 6718 in/sec 
W. “ H2 


Inter>b1ade 
Phase angle* 
0 d<!:g. 

Reduced 

Frequency, 

•^s.ref 

Frequency, 
f Hz. 

Damping, 

9 

Non-dim. 

Frequency, 

Non-dim. 

damping, 

M ® (g/2)»v 

0 

.276 

202.7 

-7.971E-2 

1.079 

-4.299E-2 

36 

.276 

203.1 

-6.950E-2 

1.081 

-3.756E-2 

72 

.284 

208.5 

-1.454E-2 

1.110 

-8.067E-3 

108 

.286 

210.2 

-2.718E-2 

1.119 

-1.520E-2 

144 

.288 

211.4 

-1.666E-2 

1.125 

-9.372E-3 

180 

.287 

211.0 

-1.339E-2 

1.123 

-7.518E-3 

-144 

.286 

210.3 

-6.401E-3 

1.119 

-3.582E-3 

-108 

.282 

207.4 

4.289E-4 

1.104 

_ .. 

2.367E-4 

-72 

.279 

205.6 

-9.847E-3 

1.094 

-5.387E-3 

-36 

.269 

197.6 

-1.625E-2 

1.052 

-8.544E-3 
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TABLE 15. ROUT LOCUS RESULTS FOR FIRST BENDING MODE, CASE 4, 1st ITERATION 


filter-blade 
Phase angle. 
a deg. 



Reduced 

Frequency, 

•^s.ref 


Vj^ref ° 


W. 


Frequency, 
f Hz. 


Damping 

9 


Non-dim. 

Frequency, 

V “ f/f 


Non-dim. 
damping, 
u ° {g/2)*v 


.258 

190.7 

-6.859E-2 

1 .065 

1 

-3.654E-2 

.259 

191.4 

-6.096E-2 

1.069 

-3.259E-2 

.265 

196.0 

-1.769E-2 

1.095 

-9.685E-3 

.267 

197.3 

-2.765E-2 

1.102 

-1.524E-2 

.268 

198.1 

-2.444E-2 

1.107 

-1 .352E-2 

.268 

198.1 

-1 .836E-2 

1.107 

-1.016E-2 

.267 

197.1 

-1.402E-2 

MOl 

-7.719E-3 

.265 

195.6 

-8.329E-3 

1.093 

-4.551E-3 

.261 j 

193.1 

-9.617E-3 

1.079 

-5.187E-3 

.254 j 

187.9 

-7.579E-3 

1.050 

1 

-3.978F.-3 
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Figure 63. SR-5 CLASSICAL FLUTTER SUMf^ARY 
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Analysis: Kaza, Kielb (Ref. 7) 

Beam modes, 2-d sub. case. aero. 



5.0 CONCLUSIONS 


1. The two-dimensional subr.onic cascade unsteady aerodynamics theory of Ref. 1 
has been 

a) modified to account for the variable sweep angles of the blades of 
the advanced turbopropellers, and 

b) Implemented In NASTRAN Level 17.7 for use In a strip theory manner on 
blade chords selected normal to a "spanwise" reference curve. 

2. The blade has been structurally and aerodynamically modelled using existing 
and modified NASTRAN finite elements capability. 

3. The In-vacuo natural frequencies and mode shapes of the blade. Including 
the differential stiffness effects due to centrifugal loads, have been 
computed. 

4. The two-dimensional strip aerodynamics and the structural modal properties 
have been formally Integrated to determine the generalized aerodynamic 
coefficients matrix for the blade modes. 

5. Modal flutter analyses conducted on two SR-5 five- and ten-bladed propellers 
have shown good correlation with NASA wind tunnel tests. 


6.0 RECOMMENDATIONS 


This is a new capability 1n NASTRAN to conduct modal flutter analysis of 
advanced turbopropellers with variable sweep blades and, although Its 
validity has been clearly demonstrated by the limited number of examples 
presented, It should be exercised extensively to determine Its potential 
and limits. 

Development and application of a three-dimensional unsteady compressible 
aerodynamics theory for multl-bladed propellers, In a manner similar to 
that of Ref. 8 , Is a desirable extension of the present work. This should 
ameliorate the unknown effects of "acoustic resonances" present In the 
two-dimensional approach. 


APPENDIX A 


GENERALIZED OSCILLATORY AERODYNAMIC COEFFICIENTS MATRIX Q 
FOR THE SWEPT BLADES OF ADVANCED TURBOPROPELLERS 


A.1 GENERAL 

The generalized modal aerodynamic force matrix [Q^^l i^or flutter analysis 
of the advanced turbopropeller with swept blades Is derived. The two-dimensional 
subsonic cascade unsteady aerodynamics program of Rao and Jones (Ref. 1 ) has 
been modified to Include the effects of blade sweep. The blade Is spanned by a 
number of non-intersecting chords selected normal to any spanwise reference curve 
such as the blade leading edge (Figure 64). The modified two-dimensional cascade 
theory Is applied on each of these chords to determine the generalized aerodynamic 
forces acting on the associated strip. The strip results are added to obtain the 
blade aerodynamic matrix. 

In order to conform to the existing computational scheme, and thus reduce 
extensive recoding, the chordwise generalized aerodynamic matrices are first 
computed for predefined, simple, chordwise aerodynamic modes (Figure 65). The 
chordwise structural modes are then Introduced via Fourier expansions to finally 
obtain the blade generalized airforce matrix, 

A.2 DEFINITIONS 

Figure 64 illustrates some of the definitions pertinent to Incorporating 
sweep effects In the 2-d cascade program. A_B_, AB and A^B^ represent three 
successive chords with points A's on the leading edge. For the chord AB, at 
any operating condition represents the absolute inflow velocity while AlT 
(= J? X lUT) is the blade (tangential) velocity. WA and AU uniquely define a plane in 
which the Inflow properties are defined. 

In the plane WAU, VA = 15 a - AlT represents the relative Inflow velocity. 

AI Is the line of Intersection between the axial plane through point A and 
the plane WAU. Angle lAV defines the relative inflow angle 6 (shown positive). 
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The angle of sweep A Is defined as the angle of Inclination of the chord BA 
with the plane MAU. A shown in Figure 64 Is positive. 

AD is the projection of AC (BA extended to C) In the plane WAU. Angle lAD 

represents the stagger angle X, and Is shown positive. 

A local coordinate system xy5 Is defined at the leading edge point A of the 
chord AB such that x is directed along AB. y is defined normal to the 'mean' sur- 
face containing the points A_, A, A^, B^, B and B_. The unit vector along y, for 
the sense of n shown In Figure 64, is given by 


! 


1 

I 


(rB^) X (AB) 

j(aX) (^1 


(AB) X (A^BJ 
i(AB) X (A^BJI, 


(Al) 


Modal translations along y and rotations about x are used in deriving the 
generalized airforce matrix (Section A. 6). For the opposite sense of rotation, 
xyz is defined to be left handed with y reversing diroction. 

The shaded area about the chord AB represents the strip of integration 
associated with AB. The length of the strip Is given by AB, the widths at the 
leading and trailing edges are respectively given by 

e. “ 2 ^ 

“t.e. = i (M*) • ^ . 

A 

where k Is the unit vector along i. 

A. 3 BLADE DEGREES OF FREEDOM 



For aerodynamic computations, consider the propeller blade to be spanned by 
S chords with G structural grid points on each chord. The blade degrees of free- 
dom (u) can then be partitioned as 


iu) =< 


“1 
“2 

"s 


(A3) 


S J 
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“11 


“52 


{U|) : } 


“sg 


“sG 


(A4) 


Ujg represents the 3 translational and 3 rotational degrees of freedom at grid 


point g on chord s. 

A. 4 NORMAL DISPLACEMENT. DOWNWASH AND PRESSURE DISTRIBUTION IN TERMS OF CHQRDWISE 


AERODYNAMIC MODES 

The normal displacement (along + y) at any point x on the chord i can be 
expressed In the local (chord) coordinate system as a linear combination of the 
chordwise aerodynamic modes: 


u®"(x, z, t) = i-{z) L*^(xLl|n^^(i)| e 


i(i)t 


(AS) 


where 4; is the semi-chord, 


- L . * 


2S, 




2IL: 


represents the chordwise aerodynamic mode shapes, and 


f '1 


; f- 

a 


'1 


‘2 

‘g 


V J 


(A7) 


is the aerodynamic modal participation vector. 
The downwash at x can be written as 


t) - + V; ( 4 - r taoA ^ >] 1 . 1 ) . 


S '3x 


3Z 


(A8) 


1 
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where the chordwise cascade relative Inflow velocity from Figure 64 Is 

V| » [(VA) cos (6-X)l cos A “ V cos A . (A9) 

The quantity within the brackets represents the streamwise cascade relative 
Inflow velocity. 

Substituting for u®” from equation (A5)^the downwash can be written as 
z, t) » V|[|,(ik + (-§y^ tan A) GL ^ 

(AlO) 

©tan AL<>^J{n^^)3e^‘^'^ , 


where 


and the reduced frequency 

k = cos A /{V cos A) = w£/V . (12) 

Corresponding to the downwash boundary conditions associated with each of 

il s2 

the elements of the modal participation vector (n } and its derivative {n )i 
the modified 2-d subsonic cascade unsteady aerodynamics program computes the 
differential pressure at P discrete points distributed along the chord s. (The 
contribution due to the variation of the perturbation velocity potential in the 
z direction is neglected.) The physical differential pressure at these P points 
can then be found from 

(A13) 




si 


3 (It) 


(All) 
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In order to properly account for the leading edge singularity In subsonic 
flow while Integrating the pressures along the chord, the modal differential 
pressures are rewritten 


- [n^pj [cJp]T , and 

- [n„J [cJjT 

pg pp gp . 


where 



?T*tr*x, 

sin ( ^ 

2iz 



sin 


P-1 

( - 



) 


(A14) 




sin 


( 


sin 


^P-l-TI’Xg 


) 



(A16) 



sin 


FT.ff. 

( 



sin ( 


P-l-ir-Xp 


za- 


) 


A. 5 GENERALIZED AIRfORCE MATRICES FOR CHORDWISE AERODYNAMIC MOOES 

The differential pressures are estimated at P points along the chord (equa- 
tion A13). The physical displacement at these P points can be written using 
equation AS as 



(A16) 


The virtual work done by the aerodynamic forces, acting on the blade strip 

snn 


associated with the chord s, through a virtual displacement 6{u } becomes 

ou ' P® • (C^’x + 

X“0 


6W® = / ®6u® 


dx 


(A1Z) 


OfJ!G«MAL’ Pf OE 13 
OP POOR OHAUTY 


r«. 4 P| v-2 ^.3 6Ln=b[*Jgf r M 0 [np^ { [cjp3 (n^h ^ [C^pf {n^2j) . 

{A18) 


where 


fuj 


ji Xi -s2 Ai 

[c' ( IT ) + <ir)J 


,s2 


[C*'( r?- ) ♦ [''(k)J 


o 


o 


5i 

Cc'^Tt) ^ V] CM-^}] 


a-" ■■ 

S S 


t$i D-u^ J + ^ f J 


Substituting for W 3 Ip equation (A18), 


6W 


' ■ ? Os ( [''ggAn*') * ) 


gg- 


99* 


whereby -i*? generalized chordwise aerodynamic force matrices are 


^%g^ “ ^^gp^ ^^gp-* * 


[A'^] » 11%: . 


with 


99* 

pS 


C*2 ..1 


[ijpj 


tl^p] = r J [HppJ 


and 


[ijpj = f w' g [Hpp] 


(A19) 


(A20) 


(A21) 


{A22J 


(A23) 


2 1 

Elements of matrices Igp and Igp have been evaluated as chordwise integrals, and 


are presented in Appendix B. 
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A. 6 TRANSFQPMATION BETWEEN AERODYNAMIC AND STRUCTURAL MODAL COORDINATES 

Equation (AS) can be used to express the normal displacement at the 6 struc- 


tural grid points on chord s as 

tu®"} = i^sC^ggKn^b . (A24) 

In terms of structural modal coordinates, can also be written as 

(u®"} » . (A2S) 

Comparison of equation (A24) with equation (A25) yields 

tn“’> “ i U,) . (A26) 

Differentiation with respect to z/i- results In 

0 



where 


r«s2-, r*i -i"^ r 3 Jn.locali 

[6glJ = [4>ggJ <t>gi J 


(A29) 


The modal displacement matrix (jjSf'jlocal derivative ^sn, local 


derived from the global modal matrix 4) as follows; 


[4?’—] = wljl] 


(A30) 


where each of the 1 modal columns of consists of the three modal 

translations at each of the G computing station grid points on the chord s. 

enl 

The transformation G Is given by 
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0 1 0 
0 1 0 

0 1 0 


[[rf] i 


[T|?f 




O 


O 


[T^9]T 
«• sgJ 




(A31) 


where and T*’® represent the coordinate transformations from the NASTRAN basic 

system to the local (chord) xyi system and the £loba1 (displacement) system, 
respectively. 

Similarly, 

r 3 iSn,locali rpinZ-i r .s2,Qlobal-i 

where each of the i modal columns of ^^2, global ^Q^s^sts of the three modal 
rotations at each of the 6 computing station grid points on the chord i. 


The transformation 6®*^^ is given by 

1 0 O' 

1 0 0 




1 


.1 0 0 . 




tTf^] 


I 






'■ sgj 


'TT^g-iT 

'■‘sG-' 


where m-j = 1 for right handed xyz, and 
= -1 for left handed xyz. 


} 


(A33) 

(A34) 


A. 7 GENERALIZED AIRFORCE MATRICES FOR CHORDWISE STRUCTURAL MOOES 

Ihtroducing relations (A27) and (A28) in equation (A21), the virtual work 
expression for chord s becomes 


1 


2 


ew' . t p-s l^l\f % [G=|] 

u- 
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whereby the generalized airforce matrices for chordwise structural modes can be 
obtained as 


[q||j " I P 5 ^ [A=p [G=;] = i Pj IQ?;], and 


2 „ 2 1_ puSlnT mill rr*Sli _ L 

's 


2 „ 2 rnSl' 


wfiJ = \ Pi "s' ‘s' ) 

“ ? P| ''s' ‘s' «|?3 

Combining and the generalized airforce matrix for chord s Is 
defined as 


(A36) 


[Q^l] ° 2 ^s.ref 


(.2 A 

's.refy 
? 

2 P|,ref ''i.ref 


i Ps ''i' ‘i' 


*^s,ref '^Sjref 


)(tq?]] * [Qf?]) . 


(A37) 



A.8 BLADE GENERALIZED AIRFORCE MATRIX 

The virtual work done by the aerodynamic forces on the blade Is the sum 
total of the virtual work done by all the chordwise aerodynamic forces on their 
respective strips. Thus 

5„bla(!e.| j,,s 
i=l 

“ Z ^s,ref ^s^ref ^ 6L^^J , (A38) 

s 


whereby the generalized airforce matrix for the swept blade of the advanced turbo- 
propeller can be written as 

[Q^jlJ ® 2 Pi, ref '^s 
s=l 


= 2 Ps,ref '^s.ref ^ ^^11^ “ 2 Ps,ref '^s.ref ^^ii^ ‘ 
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CHORDWISE WEIGHTING MATRICES 

The chordwise weighting matrices I_ and I appearing In equations (A23) 

yi" yp 

of Appendix A are evaluated. 

1. r wi'j [Hppj 



'x.l 

^x.2 

• * • ^x.(P-l) 


‘xM 

^x^2 

* • ^x^(P-l) 


'x.1.1 

^x.1.2 

’ • * ^x.l.(P-l) 

^x.2.y 

'x.2.1 

^x.2.2 

• * ‘ ^x.2.(P-l) 


Jx.(G-2)./' ^x.(G~2).l ^x.(G-2).2 ‘ ‘ ' ^x. (G-2). (P-1 )J 


The six typical elements of Igp are as follows: 


I ( 


0 

ST/ 2 


t)/ 


'2-(x/il|) 

(X/lj) 


<h) 




c) I 


xV 


4 nr' 


” Ai)'/ 


2-(x/H-) 


«(i) 


1* 2g •••• ( P“1 ) • 


= ir/2 . 
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d) 


*x^.r 


Sin d (^) 




8 

?if 


t r even, and 


= 8 32 ^ 

7S ■ (rif)> • '’‘‘‘*- 
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r°1, 2 (P-1) . 




1%) • (^) 


•(i/lj) 




M, 2. 


This Integral Is evaluated numerically, and the values are tabulated as 
follows; 


f) I 


x.r.s 


r 

^x.r.v^ 

1 

1.1333 

2 

-0. 00036 

3 

0.18796 

4 

-0.00027 

5 

0.08469 

6 

-0.00022 

7 

0.05049 

8 

-0.00019 



s 

1 s 2y • ■ • t 


1) 

^ . 


1, r=s 

0, r^s, (r+s) even. 


- r(r-f)^‘|Wsr • 


(S-2) . 
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2. [ijpl " f W’j [HppJ 


ORlCiNAL PAGE IS 
OF POOR QUALITY 


■V 

'l. 

^2. 

• • • I(p-i) 

■XnA 

'x.l 

'x.2 

• * ■ ^x.(P-l) 

I.-/- 

'l.l 

^1.2 

• • • M.(P-l) 

Z.V’ 

*2.1 

^2.2 

' * • ^2.{P-1) 


(G-2)V^ 


rl 


I 


(G-2).2 


The six typical elements of I^p are as follows: 
2 /2-(;/ll;) , - . 

’’ ^ I /-wqr '> i-q) 


I 


(G-2).(P-1) 


») >r “/ • -• 

■ ^[1 - (-If]. 


= 1.2 P-1 . 


c) 

d) 

e) 


n 

IxV^as given by typical element (a) of Igp. 


Iv .. given by typical element (b) of . 

X # r up 


Ip/-= / sin 
^ 0 


/ r«nr*x > / 


gp 

W^T ^ V n- y * 


r = 1, 2, .... G-2* 


This Integral is evaluated numerically, and the values are tabulated as 


follows; 




» 


ORIGINAL PAGE iS 
OF POOR QUALITY 


r 

V.V" 

1 

1.48283 

2 

0.89414 

3 

0.83521 

4 

0.66721 ' 

5 

0.64172 1 

6 

0.55519 1 

7 

0.54026 

8 

0.48547 



r = 1( 2f !••• G-2 I 
s = 1 , 2» . . . , P-1 . 

= 1 , r=s 

= 0 , rH 
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SYMBOLS 


A Chordwise aerodynamic matrix for aerodynamic modes 

B Blade spacing 

C Matrix of modal pressure coefficients, constants 
c Chord 

f Frequency (w/2it) 

G Transformation matrix 

g Artificial structural damping 

I Chordwise weighting matrix 

i /=T 

K Stiffness matrix 

k Reduced frequency 

H Semi chord 

M Mass matrix, cascade relative inflow Mach number 
m^ Constant 

N Number of blades on propeller 

P Load vector 

p Pressure, complex eigenvalue 

Q Aerodynamic coefficient matrix for structural modes 
T Transformation matrix 

t Time 

u Physical degrees of freedom 

V Cascade relative inflow velocity 

W Virtual work, weighting matrix 

w Strip width, downwash velocity 

7, y, z Chord local coordinate system 
6 Relative inflow angle 
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SYMBOLS (Continued) 


n Aerodynamic modal coordinates 

A Sweep angle 

X Stagger angle 

y Non>d1mens1ona1 damping (gv/2) 

V Non-dimensional frequency (^/ ^vacuum) 

^ Structural modal coordinates 

p Mass (flow) density 

a Interblade phase angle 

^ Chordwise aerodynamic modes 

<J) Structural modes 

T Differential pressure distribution In aerodynamic modes 

Q Rotational speed 

u) Circular frequency 

SUBSCRIPTS 

g Grid point on chord, chordwise aerodynamic mode 
1 Structural mode 

J2,.e. Leading edge 

p Cascade theory pressure point 

ref. Reference 

s Chord s 

t.e. Trailing edge 


SYMBOLS (Continued) 


SUPERSCRIPTS 


b Basic coordinate system 
cf Centrifugal 
d Differential 
e Elastic 


g 

global 


^Global coordinate system 


local 

n 

s 


|Local coordinate system 

Normal to chord 
Chord s 


V 
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